home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 2002 November / SGI Freeware 2002 November - Disc 3.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / viscous.ode.z / viscous.ode
Text File  |  2002-01-08  |  689b  |  27 lines
  1. # This example simulates a bead sliding on a smooth circular wire.  The
  2. # output displays the height of the bead as a function of time.  There is
  3. # viscous damping, so the bead should settle toward the equilibrium point
  4. # (i.e. ordinate = 0).
  5.  
  6. # You may run this example by doing:
  7. #
  8. #    ode < viscous.ode | graph -T X -C
  9. #
  10. # or alternatively, to get a real-time plot,
  11. #
  12. #    ode < viscous.ode | graph -T X -C -x 0 20 -y 0 2.5
  13.  
  14. a = 1  # radius of circular wire
  15. g = 10 # acceleration due to gravity
  16. w = 10 # angular velocity of circular wire
  17. b = 1  # damping coefficient
  18.  
  19. the' = vthe
  20. vthe' = (w^2)*sin(the)*cos(the) - (g/a)*sin(the) - b * vthe
  21. the = 0.1
  22. vthe = 0
  23.  
  24. print t, the 
  25.  
  26. step 0, 20
  27.